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Abstract 

Changing the spheroidal wave equations into new Schrddinger's form, the super- 
potential expanded in the series form of the parameter aare obtained in the paper. This 
general form of the super-potential makes it easy to get the ground eigenfunctions of the 
spheroidal equations. But the shape-invariance property is not retained and the corre- 
sponding recurrence relations of the form (TJ| could not be extended from the associated 
Legendre functions to the case of the spheroidal functions. 

PACs: 03.65Ge, 02.30Gp,11.30Pb 

Introduction 

Spheroidal wave functions appear in many different context in physics and mathematics 
[l]-[3]. Their differential equations are 
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They just have one more term ax 2 than the spherical wave equations (the associated Leg- 
endre' equations). This extra term makes great difference in Equations: the spherical wave 
equations belong to the case of the confluent super-geometrical equations with one regular 
and one irregular singularities, whereas the spheroidal wave equations are the confluent Heun 
equations containing two regular and one irregular singularities. The extra singularity makes 
them one of the toughest problems for researchers to treat [T]- [3J. 

The spheroidal functions are the solutions of the Sturm-Liouville eigenvalue problem of 
Equation (prj with the natual condition finite at the boundaries x = ±1. The eigenvalues 
are the allowable values Eq, E±, . . . ,E n , ... of the parameter E; and the spheroidal wave 
functions are the corresponding eigenfunctions Oo, Oi, • • • , & n , - - - |T] - [3] . 

Though the spheroidal functions are the extension of the associated Legendre-functions 
P ; m (x), they stand in vivid contrast against the associated Legendre-functions. The associated 
Legendre-functions PJ n {x) have mane elegant properties: 

1. When m = 0, the Legendre-functions P n {x) = Pn{ x ) are polynomials of the independent 
variable x; 
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2. Back to the variable 6 with x = cos 6, the recursion relation among the Legendre- 
functions could be written as: 



P n (cos6>) = n' 1 



-ncost/ — sm( 
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P n _i(cos0), n = 1,2,3, 



so, all P n (cos9) could be deduced from the ground function Pq from the recursion 
relations 

3. With the quantity m fixed, all the associated Legendre-functions P™(cos0) could be 
derived from by 



P™(cosO) = [(n + m)(n — m)\ 2 



-n cos v — sm ( 



p^ose), 



n 



m + 1, m + 2, m + 3, . . . , . 



(3) 



4. With the quantity n fixed and the transformation P™ = t ' — the recurrence rela- 

sin2 

tions between adjoint m of the associated Legendre-functions become 



{m+ 2 )cote -Te 



*„(0,m), 



* n (0,m + l) = [(n + m + l)(n - m)} 2 * 

m = 0,1,2,... ,n. (4) 
These equations @ are more useful than that of ([3|), they result in the relation 

™ = (l-, 2 ) f ^^. (5) 

The associated Legendre-functions -P™(x) are related to the Legendre's functions by ([5]), 
which is the repeated forms of the relations (HJ). The significant meaning of these relations 
are that one only needs to know the Legendre's functions to obtain the associated Legendre's 
functions. So the two important works of the associated Legendre's equations lie in finding 
the ground state function of the Legendre's equations, that is, the case of m = and the 
recurrence relations of the type of Eqs.([3]) and dH). Eqs.© could be used to get the excited 
state eigenfunctions of the Legendre's equations; while Eqs.(dD make it possible to obtain the 
associated Legendre's functions from the Legendre's functions. 

Like the spherical wave functions, the spheroidal wave functions are important functions 
applied extensively in many different branches in physics and mathematics PQ-[1], hence, 
it is natural for one to require whether or not the similar relations exist for the spheroidal 
functions. This is a long standing problem. However, it has actually been treated in the series 
papers[8]-[9]. In these papers, the supersymmetric quantum mechanics (SUSYQM) is first 
applied to study them, and new interesting results are obtained [5]-|lL)j: (l)the general form of 
the the ground eigenvalue and eigenfunction are given, which reduce to the ground eigenvalue 
m(m + 1) and the ground eigenfunction P™ when a = 0; (2) the generalized recurrence 
relations like that of © are obtained for the spheroidal functions [10]. See references [6]-[T0] 
for details. As stated before, the relations (jl]) are crucial for one to obtain the associated 
Legendre's functions from the Legendre's polynomials. Therefore, it is natural for one to 
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investigate whether the same kind recurrence relations could be extended to the case of the 
spheroidal functions. This is just what the present paper tries to do. 

As done before, the super symmetry quantum mechanics is used to deal with the prob- 
lem. The recurrence relations turn out to be the shape-invariance relations for the corre- 
sponding differential equation. Thus, the spheroidal equations are first transformed into 
the Schrodinger equations. Then the super-potential W is introduced and solved in the se- 
ries form of the parameter a. Thirdly, the ground state eigenfunction is obtained. Finally, 
the shape-invariance property is checked, and the recurrence relations are studied. In the 
following, the same steps will proceed too. 

In use of the perturbation methods in super-symmetry quantum mechanics, it is necessary 
to rewrite the differential equations in the Schrodinger form. In the previous papers, the 
original form 
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is used to study. It is obtained from Eq.flT]) by the transformation x = cos 9 and the corre- 
sponding boundary conditions become finite at 9 = 0, it. 

There are two ways to transform Eq.([S]) or ([T]) into the forms of the Schrodinger equation. 
The first transformation is to change the eigenfunction through the transformation O 

and the Schrodinger's form becomes 
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* = 0. 
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The corresponding boundary conditions now are ^(0) = ^/(tt) = 0. 

In the series papers [6]- [10], the supersymmetric quantum mechanics (SUSYQM) is first 
applied to treat Eqs.(J7J). The focus is the super-potential W, which is determined by the 
potential 

,2 1 ' 



V(0, a, m) 
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The ground function is connected with the super-potential W by 



iVexp 



Wd6 



(8) 



(9) 



Hence, whenever the super-potential is known, the ground state function is known too [6]- 
|10| . By the perturbation methods in the super-symmetry quantum mechanics, it is the super- 
potential W that is expanded in the series form of the parameter a, that is, W = J2^=o an W n . 
By this method, new interesting results are obtained [6]-|10|: the first several terms of the 
the ground eigenvalue and eigenfunction are given, which reduce to the ground eigenvalue 
m(m + 1) and the ground eigenfunction P™ of the associated Lengdre' functions. Of course, 
these results are obtained through the corresponding terms for the super-potential V7[6]-[8]: 

W = - (m + - ] cot 9, E 00 = m(m + 1). (10) 
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Later the reference [9] generalized the results and obtained the general form for the super- 
potential W = J2n=o w na n as 

Wn = y , n &n ' k ^ sin 9 cos 2fc+1 9, n = 1, 2, . . . (13) 
^ (2m + 3) n ' ' v 1 

where a n ^ could be easily determined step by step. The shape-invariance property is also 
proved in it. The form of the super-potential in (|13p is a little different in form from that 
of the references [6]. The reference [ID] directly generalized the form of the reference [B] and 
gived an alternative form of the super-potential as 



t^ n = cos sin 2 ^. (14) 

k=l 



Our motivation is to investigate whether we could extend the relations of the type (j3|) 
to the cases of the spheroidal functions. In these recurrence relations, it is the integer m 
that is different. From the point of the super-symmetry quantum mechanics, the recurrence 
relations are oriented in the properties of the shape-invariance of the super-potential and 
are the relations between different eigen-functions of correspondingly different eigen-values. 
Therefore, the integer m should stand in the position of the eigen-values, the energies. Actu- 
ally, this could not be met in Eq.© where it is the quantity E that stands as the eigen- value. 
However, there is another way to transform Eq.© to the Schrodinger form. Changing the 
independent variable 9 to the new one z = lgtan(|), the new Schrodinger 's form for the 
spheroidal wave equations is obtained as: 



d 2 e 

dz 2 



+ 



Esech 2 z + asech 2 z — asech 4 z — m 2 



9 = 0. (15) 



First, one notices that the interval (0, tt) in the original variable 9 now corresponds to the 
interval (— oo, +oo) in the new variable z. Secondly, the boundary conditions now turn out as 
O finite at z — > ±oo. The most important thing is that Eq (|15|) makes the term containing the 
original eigenvalue E no longer in the position of the eigenvalue. Instead, it is the quantity 
—m 2 now that is the eigenvalue of Eq. (|15|) . This is just what we want to obtain. It will 
beneficial for one to further investigate the kind relations 0] for the spheroidal functions. 
In the above Eq. ([T5|) . the potential is 



V(z, a) = — -Esech 2 ,? + asech 2 z — asech 4 z . (16) 

When a = 0, Eq. (|15p is just the form from the associated Legendre equations ©; its 
ground energy is — m 2 , and the nodeless ground eigenfunction is the associated Legendre 
function with the original eigenvalue E taking m(m + 1). Therefore, when a ^ 0, 
the ground energy is also — m 2 and the ground eigenfunction is Q(z,m, E), which must be 
nodeless. Actually the original eigenvalue E could be determined in the following by requiring 
the eigenfunction finite at the infinities. Its value could also be obtained from the previous 
papers [9], [TO] . 

The super-potential W is determined from the potential V(z, a) by 

W 2 -^- = V(z,a)+m 2 (17) 
dz 
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where subtracting the ground eigenvalue — m 2 to make the potential factorable. This equation 
is the same difficult to treat as the Schrodinger original Eq. (fT5j) . hence, the perturbation 
method is used to solve it. First, the super-potential W is expanded as the series of the 
parameter a, that is, 



W = W + aWi + a 2 W 2 + a 3 W 3 + ... = £ a n W n . 



(18) 



ra=0 



W 2 - W = Wl — Wq + q (2WWi - W[) + a 2 (2W W 2 + Wl - w£\ 

+ a 3 (2W W 3 + 2W 1 W 2 - W 3 ) + a 4 (2W W A + 2W X W 3 + W% - W$\ + . . . . (19) 
Secondly, the original eigenvalue E must also be written as 



a 



(20) 



n=0 



where there are three lower indices in the parameter -Eo,n;m with the index refereing to the 
original ground state energy, and the other index n meaning its nth term expanded in the 
series in parameter a and the last m indicating the parameter —m 2 in Eq. (|15p . One can 
write the perturbation equation as 



W 2 - W' 



V(z, a, m) + m A 



-E'sech 2 ^ + asech 2 z — asech 4 ; 



+ m 



a n .E'onSech 2 z + asech 2 z — asech 4 



.n=0 



+ m 



Comparing Equations (fT7|l . (fl~9j) . and (f22|) . one could get 

Wq -Wq = -E 0t0]m sech 2 z + m 2 

2W W! - W{ = sech 4 z - (£ ,i;m + l)sech 2 z 

2W W 2 + Wl -W 2 = 

2W W 3 + 2W X W 2 - W 3 = 



-£j ,2;msech z 

"£fl.?!:mSech 2 2I 



(21) 
(22) 



(23) 
(24) 
(25) 
(26) 



n-l 

2W W n + Y / W k W n - k -W^ = -E Qtn . m sech 2 z 
fc=i 

From Eq. ([23|) . we get 

Wo = mtanhz, -Eo,0;m = m(m + 1). 
Then, we can write the other equations more concisely 

W[ — 2m tanh zW\ = (E ^- m + l)sech 2 z — sech 4 z 
W 2 - 2m tanh zW 2 



W 3 - 2m tanh zW 3 



E 0t2]m sedi 2 z + W 2 

£ ,3;mSech 2 Z + 21^1^2 



(27) 

(28) 

(29) 
(30) 
(31) 



WL - 2mtanhzW n 



n-l 
fc=l 



(32) 
(33) 



5 



After obtaining the zero term Wo for the super-potential W, the first order W\ can be gotten 
as 

Wi = A x cosh 2m z (34) 

with 

dA x 



de 



sech 2m z 



(-E-o,i;m + l)sech 2 z — sech 4 ,z 
(■Eb,iim + l)sech 2m+2 z - sech 2m+4 



dz 



(35) 
(36) 



In order to simplify the calculation, we just write some useful formula [19 

Q(m,z) = 



sech 2m zdz 



sinhz 



2m - 1 

This formula could be written in concise form 



sech 2 - 1 * + V 2 fc (m-l)(m-2)...(m-fc)sech"- 1 
j^J (2m-3)(2m-5)...(2ro-2Jfe-l) 
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Hence, 
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(2m - \){2m - 3)(2m - 5) ... (2m - 2/c - 1) 



Q(m + n,z) = J S ech^ m+n ^ z dz 



sinh z 
m + n 



m+n— 1 



^ /(m + n, /c)sech 
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Qi(m + n, z) + 2(m + l)I(m + n, n — 2)Q(m + 1, z) 



where 



Qx(m + n,z) 



sinhz 
m + n 
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See the appendix 1 for details. Hence, the quantity A\ is 

A~i = (E 0i i ;m + l)Q(m + 1, *) - Q(m + 2, z) 

= (E 0>1;m + l - 2(m + + 2, 0)) Q(m + 1, 2) - Qi(m + 2, z) 

= aiQ(m + 1,2;) — Qi(m + 2, z) 

where the introducing a\ = £?o,i;m + 1 — 2(m + l)/(m + 2, 0) is used. So 



Wi = ^icosh 2m z 



axQ(m + l,z) cosh 2m z - Q x (m + 2, z) cosh 2m 2; 



(37) 



(38) 



(39) 
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The quantity £?Q,i ;m needs to be determined by the boundary conditions that the is finite 
at ±00. This in turn require to calculate the ground eigenfunction upon to the first order by 



iVexp 
iVexp 

Nsech. m z exp 



Wdz 



J W dz - a J W x dz 
a [ W x dz 



* exp 0(a 2 ) 
* exp 0(a 2 ). 



(45) 
(46) 
(47) 



Whenever the eigenfunction is obtained, the boundary conditions finite ^±00 would choose 
the proper £?o,i;m- it is easy to compute 



J Widz = J axQ{m + 1, z) cosh 2m zdz -J Qi(m + 2, z) cosh 2 " 1 zdz 

By the relations (f3"5j) . (pfT)) . (f4"2"j) it is easy to obtain the eigenfunction 

I(m + l,k) cosh 2fc z 



(48) 
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A r sech m zexp a a\ 

V fro 2fc (- + !) 
I{m + 2,0)sech 2 z 



: CXp 



a- 



2m + 4 



exp 0(a 



(49) 



Unless the coefficient a\ of the term J2T=o ^^^{m+i^ — ~ becomes zero, the eigenfunction 
will blow up at infinity whenever a%a < 0,. So, £?o,l;m is determined by a\ = £fa,l;m + 1 — 
2(m + 1) J(m + n, 0) = 0, the results are 



^0,l;m 

6o 



2m + 3 



-Qi{m + 2, 2) cosh 2 " 1 z = _ smn ^ j( m + 2) o)sech 3 z 



sinh zsech z 
2m + 3 



./Vsech" 1 ^ exp 



sech 2 z 
4m + 6 



a- 



m + 2 



■ exp 0(a 2 ). 



(50) 

(51) 
(52) 



With the first order term of the super-potential W\ , one could compute the second term 
W2 by the same process. 



Wo = Ao cosh 2 " 1 z 



(53) 



with 



A 2 = J sech 2m z (^ ,i;msech 2 z + W?) dz 
= E 0) i- m Q(m + l,z) + J sech 2m 2i 



sinh 2 zsech 6 z 



dz 



%;mQ(m+ l,z) + 



(2m + 3) 2 



(2m + 3) 2 
(Q(m + 2, z) - Q(m + 3, z)) 
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(2m + 3) 2 

[Q 1 (m + 2,z) -Qi(m + 3,z)] 



Q(m + 1, z) 



(54) 



(2m + 3) 2 

Similarly, the coefficient 02 = £?o,i;m + (2m+3)' 2 (H m + 2,0)— I(m + 3, 1)) of the term Q(m + 
1, z) must zero. The results are elegant: 

2m + 2 



0,2;m 



w 2 



"(2m + 3) 3 (2m + 5) 

sinh zsech 3 z sinh zsech 5 z 

(2m + 3) 3 (2m + 5) ~ (2m + 3) 2 (2m + 5) 



(55) 
(56) 



We can continue to calculate W3, W4 



, ...,w n , 



here we just use induction to prove 



that 



W n = sinh z 53 o n) fcsech 



2/c+l. 



(57) 



fe=i 



with a n .fc could be determined by ajj,i <n,j < n, which appear in Wi, i < n. Obviously, 
Wi in Eq. ([5T]) satisfies the requirement (f57|) in the case n = 1. Suppose that Wk, k < n 
meets the requirement, then 

W n = A n cosh 2m 



with A n is determined by 

An 



sech 2m z 



n-l 



,n;m 

sech 2 z + ^ W k W n - k 



k=l 



dz. 



(58) 
(59) 



The first thing is to simplify the following term: 
n— 1 n— 1 n— k 

53 W k W n - k = 53 £ £ a k,i a n-k,j sinh2 ^sech 2(4+j)+2 z 
fc=i fc=i i=i i=i 

n— 1 n p— 1 

= 53 53 53 a k, P ~ja n -k,j sinh 2 zsech 2p+2 z 

k=lp=2j=l 
n— 1 n p— 1 

= 53 53 53 a k, P -ja n -k,j sech 2p z - sech 2p+2 z 

fc=lp=2 j=l 
n— 1 n p— 1 

— 5 ■ 5 < ^ y [^fe,p— jQ>n— k,j Q'k,p—l—jQ"n—k,j\ Sech P Z 
fc=l p=3 j=l 

n—l 2 n—1 n 

+ 53 53 afc,2-ja n _ fc jsech 4 z - 53 53 afc,n-j«n-fcjsech 2n+2 z 

fc=l j=l k=l 3=1 

n+1 

= Wech 2p z (60) 

p=2 

n— 1 p— 1 

Vp = 53 53 i a k,p-j a n-k,j — 0>k,p-l-j a n-k,j] , P = 2, 3, . . . , n + 1 (61) 

fc=l i=l 
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with the conditions i < j, orj < 1 then 



Oij = 



(62) 



n+l 



A„ = E 0> n;mQ(m + I, z) + ^ b n,pQ(m + P, z) 

p=2 

n+1 m + 1 

+ 2 \ 6 n „ ■ — I(m + p,p-2) 

m + p 

p=2 l 

n+l 

+ b n ,pQi(m + p, z) 

p=2 



Q(m + 1, z) 



(63) 



In order to make the eigenfunction finite at infinity, the quantity £o,m,n+i is selected by 



that is, 



n+l 



%;m + 2^ym + l)/( m + p,p-2) = 0, 

p=2 



hence, 



n+l 

E ,n;m = ~ 2 Y1 b n,p{ m + l)I{m + p,p - 2). 
p=2 



n+l 

W n = A n cosh 2 " 1 z = ^2 b n, P Qi (m + p, z) cosh 2m z 

p=2 

sinh. z ? 

= J2 b n, P ; — Y] I(m + p,p- 1 - j>ech 2j+1 z 

p=2 r 3= i 

n n+l T ( . ■• -\ 

■ , \r- \r- , I{m + p,p-l-j) 2 j+i 
= svahz } j > y »n,p sech JT - 



j=ip=j+i 



= sinh z a n jsech 2 - 7+1 z 



(64) 



(65) 



- A I(m+p,p - 1 - j) 

p=j+l 



m+p 



(66) 



we see that a n j are determined by which are completely defined by aij,i < n,j < n. 
this completes our proof. Therefore, the super-potential W could be written as 



W = W Q + ^W n a n = VF" + ^sinhz^a nii sech 2j+1 za n (67) 

n=l n=l j=l 



The ground eigenfunction becomes 



0o = N exp 
= N exp 

= iVexp 



Wdz 



/OO p 
W dz -^2° n / W n dz 

n=l 

/„ oo n 
mt&nhzdz — / a n a n j sinhzsech 2 -' +1 z(iz 



n=l j=l 



i ^1 2 i 

n=lj=l ■* 



= A r sech m 2:exp 

Back to the independent variable 6, sechz = sin#, the above ground eigenfunction becomes 



oo n n 



i ^ 2 i 

n=l 7=1 J 



G = iVsin m #exp 
The original ground eigen-energy is 

oo 

E . m = m(m + 1) + XI Qn ^o,n 
with Eo n - m is determined by Eq.([6 



sm 



2 i, 



(69) 



(70) 



n=l 



Comparison with the former results 

The elegant forms ([69|) . (|70p .(|64| ) for the ground eigenfunction and eigenvalue are the same as 
that in the reference [10J. There appeared two forms for the super-potential: W(z) of the 
form (|67p and W{6) of (|14p in the reference [10] . Whatever forms may the super-potentials 
be, they should give the same eigenfunctions for the spheroidal equations. Thus, it results 
the relation between them as 



e o (0) = J w(z)f e de. 



By z = lgtan |, it is easy to get 



W(0) = W(z)^ = -^-W(z) 



dO sin 6 

Writing W(z) back as the function of the original independent variable 9, one gets 



71=1 



W{6) = W Q {6) + Y,W n {e)a 
W (6) = -mcotO 
W n (6) = W n ~- 



cos v > a n ^ sm 
fc=l 



,2fc-l 



(71) 

(72) 

(73) 
(74) 
(75) 



10 



Comparing the results with that ()14|) of the results of the reference [TO] , they are the same 
except for the Wq(9). The difference between them is originated from the eigenfunction's 
relation 

e = (76) 

sin 2 9 

= / W{9)9. (77) 

where W{9) are the forms of (|14p . The advantage of the method is that it can be easily 
extended to the spin- weighted spheroidal equations. It will be our further study. 



Obviously non-shape-invariance property of the spheroidal func- 
tions 

In order to check whether the spheroidal wave equations have the shape-invariance property, 
the super-potential W is rewritten as the following form: 

oo n 

W(A n j, z) = Aofim tanh z + sinh z ^4 nj a n jsech 2j,+1 za n , (78) 

71=1 j = l 



With the definition 



W n (Aij,z) = sinh z A n ja ntj sech 2j+1 z (79) 

3=1 

Then, V^ z (A n j) are defined as 

oo 

V ± (A nij ,z) = W 2 (A nij , z )tW' = J2 a n V±(A itj ,z). (80) 

n=0 

We will check whether or not (Aij , z) are related with by the relations 

V+(A iJ ,z)=V-(B iJ ,z) (81) 

step by step. 

First, we write the special cases for V^, n = 0, 1. When n = 0, 

Vq = W 2 (y4 ,o, z) ~ Wo(Ao j0 , z) = m 2 Al fi tanh 2 z - A , msech 2 z (82) 

= m 2 A% - (m 2 A% + A 0>0 m)sech 2 z (83) 

y+ = Wjj(Ao,o, z) + W$(A 0fi , z) = m 2 Al fi tanh 2 z + A 0fi msech 2 z (84) 

= m 2 Aq — {m 2 A^ q — Ao t om)sech. 2 z (85) 

= V -(B 0>0 , z) = m 2 B 2 fi - {m 2 B 2 fl + £ , m)sech 2 z + Rq(A 0<q ) (86) 

i?oOVo) = m 2 A 2 - m 2 B 2 (87) 



Therefore 



£ ,o = A),o - -, #o(A),o) = 2771,40,0 - 1. (88) 
m 
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This results is exact when a = , and it just shows that the associated Legendre equations 
have the shape-invariance properties. It is easy to get the recurrence relations (|3D for the 
associated Legendre functions from the the shape-invariance properties. What interests us 
most is whether this property could extend to the spheroidal functions, the case 
When n = 1, 



Vr(Ac, 0t Ai tl ,z) = 2Wo(A 0)0 )W 1 (A 1)1 ,z)-W[(A 1)1 ,z) 

2mA nA 1 A . , 2 , 4 A 



sinh 2 zsech 4 z + 



2m + 3 2m + 3 



sinh zsech 3 z 



1 1 



(2mA 0>0 + 3) Al , t + (2-2^,0)^1,1^ (gg) 



2m + 3 2m + 3 

F+^co,^,!^) = 2W (A 0>0 )W 1 (A 1;1 ,z)+W{(A 1 , 1 ,z) 



2mA 0>0 A 11 2 4 A h i 
sinh zsecn z + 



sinh zsech 3 z 



2m + 3 2m + 3 

(2^,0-3)^ 4 (-2-2mA),oMi,i sech ^ 

2m + 3 2m + 3 v ; 



If we require the first term with the shape-invariance property, that is, 
V+(A 0>0 ,A ltl ,z) = V^iBo^B^z) 



(2mB 0fi + 3)B 1:1 4 (2 - 2m%)B M a 

= ! ^sech z H ! ^sech z, (91 

2m + 3 2m + 3 v ' 

there exist two conditions for the quantity £?x,i to meet: 

(2m J B 0i0 + 3)5i,i = (2m^ ,o - 3)^i,i (92) 
(2 - 2m5 ,o)£i,i = (-2 - 2rai ,o)ii,i. (93) 

It is easy to see that these two conditions are not compatible, one is led to the conclusion the 
shape-invariance property is not hold for the spheroidal functions, first, we dot not admit 
them and investigate carefully. The shape-invariance property of the zero-term (|88|) shows 
the spheroidal functions have the same property as that of the associated Legendre functions. 
The new eigenvalue of the Eq. (|15|) is the quantity — m 2 , whereas the original eigenvalue E 
now is contained in the expression i£sech 2 z. Therefore, when the new eigenvalue changes from 
— m 2 to — (m — l) 2 , the original eigenvalue E can not remains unchanged, this is the great 
difference of the spheroidal functions from the the associated Legendre functions where the 
original eigenvalue remains the same as the new eigenvalue changes from — m 2 to — (m — l) 2 . 
So the potential have no shape-invariance property. 
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Appendixl: Simplification of the quantity Q(m + n, z) 



r sinli 7 m+n ~ 1 
Q(m + n,z) = sech 2( - m+n ^ z dz = - ]T J(m + n, k)sech 2 ^ m+n ^ 2k - 1 z 

= Qi(m + n, z) + Q 2 (m + n, z) (94) 

where the two parts Qi(m + n,z), Q2(m + n, z) are: 

• u n-2 

Qi(m + ra,z) = ^^J2 I (m + n,k)sech^ m+n ^- 2k - 1 z (95) 

^ n k=0 
■ i ra+n— 1 

Q 2 (m + n,z) = V /(m + n,A;)sech 2 ( m+n )- 2 ^ 1 z. (96) 

m + n , ^ , 

k=n— 1 

In the above equation for the quantity Qi, the change of j = n — 1 — k simplifies it as 

QAm + n,z) = — V 7(m + n,ra- 1 - j)sech 2m+2j+1 z. (97) 
m + n r— J 

The second part Q 2 (m + n, z)is connected with Q(m, z). Here is the proof. First changing k 
in the above equation for the quantity Q2 to p = k — (n — 1), it is easy to obtain 

Q 2 (m + n,z) = ^^^/(m + n,p+ (n- l))sech 2m - 2p+1 z 
m + 

= ^-^^/(m+l + (n-l),p+(n-l))sech 2m - 2j,+1 z (98) 
Tn -\- n 

here the relations between I(m + l,k + /) and I(m,k) are 

7(m + l,k + l) 

2 k+l (m + l)(m + I — 1) . . . m(m — l)(m — 2) . . . (m — k) 
(2m + 2/ - l)(2m + 21 - 3) ... (2m - l)(2m - 3)(2m - 5) ... (2m - 2k - 1) 
2 * 2*- 1 (m + /)(m + / - 1) . . . (m + 1) 2 fc m(m - l)(m — 2) ... (m — A;) 
(2m + 21 - l)(2m + 2/ - 3) . . . (2m + 1) (2m - l)(2m - 3) . . . (2m - 2k - 1) 
= 2I(m + l,l-l)I(m,k). (99) 

So, 

7(m+ 1 + (n - (n - 1)) = 2/(m + 1 + (n - l),n - 2)J(m + l,p) (100) 

The quantity Q2(m + n, z) becomes 

Q 2 (m + n,z) = ^^^2/(m + l + (n-l),(n-2))/(m + l,p)sech 2m - 2p+1 z(101) 

Tft ~\~ TL „ 

■ , m 

sinn z 



2J(m + 1 + (n - 1), (n - 2)) ^ I(m + l,p)sech 2m - 2p+1 z 

m + n p =o 

2I(m + n,n-2)(m + l)Q(m + l,z). (102) 
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Therefore, 

Q(m + n,z) = Qi(m + n, z) + 2(m + l)I(m + n,n- 2)Q(m + 1, z). (103) 
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